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Abstract 

We show that the inclusion of simple anisotropic pressures stops the 
isotropic Friedmann universe being a stable attractor as an initial or final 
singularity is approached when pressures can exceed the energy density. 
This shows that the situation with isotropic pressures, studied earlier in 
the context of cyclic and ekpyrotic cosmologies, is not generic, and Kasner- 
like behaviour occurs when simple pressure anisotropies are present. We 
find all the asymptotic behaviours and determine the dynamics when the 
anisotropic principal pressures are proportional to the density. We expect 
distortions and anisotropies to be significantly amplified through a simple 
cosmological bounce in cyclic or ekpyrotic cosmologies when ultra-stiff 
pressures are present. 

PACs nos 98.80. Jk, 04.40.Nr, 04.50.+h, f i.25.Mj 

1 Introduction 

There has been strong interest in the cosmological consequences of admitting 
a 'ultra-stiff' fluid, whose isotropic pressure exceeds its energy density in the 
early stages of the universe. This situation could occur in a number of scenarios 
created by attempts to develop non-singular descriptions of spacetime which 
are applicable at arbitrarily early cosmological epochs. The simplest example is 
that provided by a scalar field, (/), with a negative potential energy, V{4>) < 0, 
in a homogeneous and isotropic universe, so that the pressure-density ratio is 



- = 7- 1 = 



and for < we can have p/ p > 1 in a regime where \V\ > ^(p . This scenario 
can arise in the so called ekpyrotic [1 , or cyclic, universe scenarios described by 
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Gasperini and Veneziano [5] and Erickson et al [3] , but can also occur in other 
theories containing effective scalar fields, for example in some compactifications 
in higher-dimensional theories [HIS]. As a result, several authors have inves- 
tigated the cosmological consequences of the presence of a simple, ultra-stiff, 
7 > 2 perfect fluid. Erickson et al [3] and Lidsey |^ showed that an ultra-stiff 
perfect fluid renders the isotropic Friedmann universe stable as an initial or fi- 
nal singularity is approached in general relativistic cosmologies. This is quite 
different to the situation when 7 < 2, where isotropic expansion is unstable on 
approach to initial and (if present) final singularities. Theorems were proved to 
establish this stability, a form of cosmic no-hair theorem, for homogeneous and 
anisotropic cosmologies, where the problem reduces to the analysis of ordinary 
differential equations [7]. The results are intuitively obvious for perfect fluids: 
there is no form of curvature or expansion anisotropy can diverge faster than 
as the mean scale factor a{t) — > 0, but a fluid with p — — l)p > p the 
density will diverge isotropically as p oc a~^'' in the same limit and dominate 
the anisotropic stresses in the a — >■ limit when 7 > 2. If a real fluid with 
rotational motion is present then the vortical energy density, fl^, goes to zero 
in the same limit, since by conservation of angular momentum [8] it evolves as 

cx a2(37-4) ^ as a ^ for 7 > 2 and Q^/p cx a^^-io ^ for 7 > 10/9. 

This stability of the isotropy and homogeneity of the Friedmann solution 
as a is a very important ingredient for cyclic cosmologies because it per- 
mits a regular transition from collapse to expansion from cycle to cycle. If 
the isotropic expansion were unstable as a — >■ then huge irregularities and 
anisotropics would accumulate and the successive cycles would be very different 
and increasingly anisotropic. In particular, attempts to follow small inhomoge- 
neous perturbations through the bounce would fail ^ . 

In this paper we will show that the studies of the likely evolution of cosmo- 
logical models as a — ^ when pressures can exceed energy density have only 
considered a restricted situation in which the ultra-stiff pressures are isotropic. 
When anisotropic pressures are also present the stability of the Friedmann solu- 
tion fails as the initial and final singularity is approached and the cyclic scenario 
loses this appealing feature. 

2 Anisotropic pressures 

Earlier analyses of the p > p cosmological problem have restricted their attention 
to the situation where the pressure is isotropic [31 However, if anisotropic 
pressures are present we should also expect the principal pressure components 
{pi,P2,P3) to be able to range over values that exceed p. This corresponds to 
a violation of the dominant energy condition, (which requires T^" > |T''''| for 
each a, 6, [10)). Physically, anisotropic pressures are to accompany anisotropic 
expansion at very high energies because asymptotically-free interactions become 
collisionless when T > lO^^GeT^ and gravitons will be collisionless below the 
Planck energy scale (T < lO^^GeV^). Thus, if there is any expansion anisotropy, 
these collisionless particles will redshift (or even blueshift) at different rates in 
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different directions and create significant anisotropic stresses with some (or all) 
experiencing pt > p. The evolution as the universe expands is complicated and 
has been studied in detail for the p < p situation, especially for the radiation 
and dust-dominated eras when p/p = 1/3 or 1, for evolution away from the 
singularity to the present [TTJ[T2]. Here, we will be interested in the behaviour 
of the cosmology in the limit of approach to the initial singularity with isotropic 
and anisotropic fluids both present. In order to establish a failure of the cosmic 
no-hair approach to the Friedmann solution it suffices to consider only the simple 
Bianchi type I universe. More complicated Bianchi type universes will provide 
further scope for anisotropic stresses to dominate but we will not include them 
here. 



3 Bianchi I with ultra-stiff perfect fluid 

Throughout this paper, we make use of orthonormal frame formalism developed 
by Ellis et al jl3| . As a simple introduction we will consider an anisotropic 
Bianchi type I universe containing a perfect fluid with ultra-stiff equation of 
state 

T^b^ P{uaub + i-f-l)gab}, (1) 

where 7 > 2 is a constant, as above. We assume the fluid 4-velocity vector u is 
orthogonal to the homogeneous space-like hypersurfaces. We use an orthonor- 
mal tetrad so that Einstein equations and Bianchi identities are put into the 
following form: 

(Tap = — SiJo-Q^ + 2e'"^^^cr^)^r2y, 

p = - a^, 

p = -3jHp, 



where 



a^^la^pa^'^. (2) 

Here, H is Hubble parameter (mean expansion rate), aap is the traceless shear 
expansion tensor, and fla is angular velocity of the tetrad frame with respect 
to a Fermi-propagated frame. We define proper time t by 

d 

and denote the time derivative with respect to t by overdot. Greek indices are 
used for space components of the tetrad. We can use the freedom to choose the 
space tetrad to diagonalise aap- The shear evolution equations then imply 

Qa = 0. 
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As is the case for the FLRW universe, these equations are not indepen- 
dent. To single out the dynamical degrees of freedom, we introduce expansion- 
normalised variables: 



Since we are interested in the evolution of the universe during its expansion, we 
also define a dimensionless time r by 

dr _ „ 



and a mean scale factor by 



We see that 



Uh. 



and therefore the initial singularity {I = 0) corresponds to r — >■ — oo. Denoting 
T derivatives by primes, we have 

n = 1 - s2 = 1 - ls„;3S'5«, 

q = 2E2 + i(37-2)Q, 
with an auxiliary equation 

n' = {2q-3j + 2)Q, 

and a decoupled equation 

H' = -{l + q)H. 
These equations can be readily integrated to give 

^ = ^ + ^tanh(^-^(7-2)r^ , 

J:^^ = ^ei(--2)^cosh-^(-^(7-2)r 

H = Foe-5(T+2^^cosh5 ('-|(7-2)Ty 
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where Epa^ are integration constants satisfying 

1 



So". - 0, 



^^Oq/sSq" — 1, 

while the constant Hq can be absorbed into the definition of time. Since 7 > 2, 
we see that E << 17 when the universe is "small", near any initial singularity 
at 1 = 0. 

If we take the general Bianchi I metric 

ds^ = -dt^ + a(tfdx^ + b{tfdy^ + c{tfdz^, (3) 
the corresponding kinematic variables are given by 



H 


1 
3 




CTll 


a 
a 


-H, 




i) 
b 


-H, 




c 
c 


-H, 




= 


(a 




= 0. 





The results for expansion normalised variables are translated into metric lan- 
guage as 



a oc 



(l-Soil)r (^-^ ^ ^i+g-3(7-2)r^"^^ 

b cx e(i'^«^^)^ (1 + Vl + e-3(7-2)r^ 
)^ (1 + Vl + e-3('r-2)^) 



3{t-2) 



3(7-2) 



When T — >■ — CX), the asymptotic behaviour is 

a, fe, c ^ oc Z — > 

and the universe becomes isotropic. On the other hand, when r is large, we 
have 

(a,&.c) - (e(l-5:oii)r^g(l-So22)r^g(l-So33)r<) ^ ;l-So22 ^ ^ I-E033 ) _ 

This is the Kasner vacuum metric with Kasner indices 

We conclude that the universe become fluid-dominated and isotropic on ap- 
proach to the singularity. This is in accord with the results of Erickson et al [3] 
and Lidsey [6]. 
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4 Anisotropic ultra-stiff Bianchi I universes 



We will now investigate how the above result is modified by the presence of 
an ultra-stiff anisotropic fluid in addition to an ultra-stiff isotropic fluid. As a 
simple model we take the total energy-momentum tensor to be 

Tab = T^b + T^ (4) 

Tab = P{Ua'U.b + {l*~l)9ab + 'Pab}, (5) 

where 7^ > 2 and Vab is a constant traceless (and has to be symmetric to be 
consistent with the metric) tensor describing the ratio between energy density 
and anisotropic stresses (note that momentum flow is not allowed in Bianchi I). 
This form for the anisotropic stress tensor is characteristic of electromagnetic 
fields and many other anisotropic stresses in the early universe (see ref [111 114) 
for a fuller discussion) although cosmological Yang-Mills fields require a different 
form of anisotropic stress tensor [TS] . The governing equations are modified and 
become 

H = -ij2_2^2_ 1(3^_2)^_ 1(3^^_2)m, (6) 
o'q/3 = -SiJo-Q/j -I- 2e^''^^f7^)^rj,, + /i'Pq/j, (7) 
p + fi = iH^-a"^, (8) 
p = -3jHp, (9) 
A = -Sj^Hp-a^pV^'^p, (10) 
Again, diagonalising shear tensor aa/s gives 

(0-33 - 0-22)^^1 + ^"^23 etc, 

as the off-diagonal parts of ([7]) are just defining relations for fla and arc non- 
dynamical. For later convenience, we introduce 

'7+ = ^(0-22 +cr33), (11) 

Cr- = ^(CT22-Cr33), (12) 

where the pre-factors are chosen so that 

These variables are sufficient to describe shear expansion tensor because they 
are traceless. The sub-case with cr_ = is the axisymmetric Bianchi type I 
metric. The remaining equations are then rewritten as 

H = -H2_2(a^ + a2_)-i(37-2)p-i(37.-2)M, (13) 

&± = -3H(T±+V±p, (14) 

p = -3-f^Hp-6{V+(7++V-(J^)p, (15) 

p = -SjHp, (16) 

= p + p + Sial + al), (17) 



6 



with the definitions 

V+ = i(7'22+7'33), (18) 

V- ^ ^{V22~Vz^). (19) 

describing the pressure anisotropies. As the previous section, we introduce an 
expansion-normahsed variable for the density: 



The equations are cast into the set 

E'i = (g - 2)E± + ST'iZ, (20) 

Vl' ^ {2q~3-/ + 2)n, (21) 

Z' = (2g-37, + 2-67'+E+-6-p_E_)Z, (22) 

q = 2(I]2^ + I]2_) + i(37-2)r! + i(37,-2)Z, (23) 

1 = n + Z + + (24) 



These equations no longer solve exactly but the qualitative and asymptotic 
behaviour of the system can be determined. 



5 The dynamical system 

The system described by equations pop to p4| is three-dimensional and com- 



pact because of the constraint (p4|) . We denote the entire system by -6^(1)- 
First, we have to find lower-dimensional invariant subsets. It turns out there 
are three which play important roles. They are defined as follows: 



Bianchi I 



Anisotropic fluid 



B{1) ^{Z^O} 



A{i) = {n = o} 

Aligned to stress 'vector' [V+^V-) 



All of these are two-dimensional subsets; B{\) and A{\) lie in the boundary 
of the whole space. If we take Z and S± as independent variables, i?A(I) looks 
like an inverted bowl with parabolic surface and 11(1) is a vertical slice passing 
through the axis of symmetry. The subset B{1) is the base of the bowl and was 
already discussed in section 3. 
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At this point, it is instructive to introduce polar coordinate for shear vari- 
ables defined by 

E+ = T,cos(j), 
S_ = Ssin(/). 

Accordingly, we parametrise the shear tensor by 

V+ = ^ — rcos^, 

V- = ~ ^ rsing. 

In these coordinates (j) measures the angle around the Kasner circle and (p — 8 
(and 4> = O + n) corresponds to n(J). The governing equations are then further 
simplified to 

(25) 

(26) 
(27) 

9)] Z, (28) 
(29) 
(30) 

The second step is to find the equilibrium points of the system and determine 
their stability properties. Time derivatives of all the normalised variables vanish 
at an equilibrium point and it corresponds to a self-similar solution. They are 
most conveniently characterised by defining the parameter 

which measures the stiffness of the anisotropic fluid compared to the isotropic 
one. 

We can now list all the equilibrium points, their metric interpretation, and 
their eigenvalues: 

FL equilibrium point defined by: 

= 1, Z = T, = 

Eigenvalues 

Ai = -3(7. - 2)a, A2=A3 = ^(7.-2)(l-a) 
We denote this by 





= (g-2)S-^(7.-2: 


)rZcos{(l)-e), 




3Z 

" ^(7* - 2)rsin(^- 


0), 


n' 


= (2q-37 + 2)f2, 




z' 


= [2g - 37, + 2 + 3(7, 


- 2)rT,cos{(t)- 


q 


= 2S2 + 1(37-2)17 + 


^(37* -2)^, 


1 


= n + z + E"^. 





8 



Kasner equilibrium points defined by: 



S = l, Q = Z = 0, <f> = arbitrary constant 

Eigenvalues 

Ai=0, A2 = -3(7, -2)(l-a), A3 = -8(7, - 2) (1 - r cos(0 - 0)) 

This set of equilibrium points, called the Kasner circle, will be denoted by 
/C. 

Anisotropic 1-fluid equilibrium point defined by: 

E = r, (1)^9., f7 = 0, Z = 1 - 

Eigenvalues 

Ai = 3(7,-2)(a-r2) 

A2 = A3 = ^(7.-2)(l-r2) 

We denote this by Ai . It becomes unphysical unless < r < 1 because of 
the constraint (|24|) and the need for positivity of S and Z. For r = it is 
identical to J- and merges into the point on the Kasner circle with (j) — 9 
when r = 1. To work out the scale factors in the same fashion as we did 
in section 3, we have to restrict ourselves to the case Vap — 0,{o- /3) 
because the metric ([3]) does not admit rotation or off-diagonal shear. The 
scale factors are 

a (X 6 oc c oc t^' 



with 



1 A 2Vu 



1 + 9 V 7* - 

2 j^- 2 2 

This has the same form as the Kasner solutions but instead of the usual 
summation relation, the exponents satisfy 

P1+P2+P3 = (32) 
1 + 9 

pI+pI+pI = ^^^(l + 2r2). (33) 
Anisotropic 2-fluid equilibrium point defined by: 
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Eigenvalues 
Ai 



r(7,-2) 



4(7,-2) 



l-a- W(l -a) 1 -9a+8— 



A3 = -(7,-2)(l-a) 

We denote this critical point by A2- For it to be physical we require 

< a < < 1 (34) 



or 



< a < 1 < r. (35) 
The metric is the same form as Ai with the Kasner exponents given by 



Pi = 

P2 = 

P3 = 

Pi + P2 + P3 = 



2 2"^ 

Pi + P2 + P3 



2 

3^ 
2 

3^ 
2 

3^ 
2 

7 

4 
37^ 



1 — a 



2Vu 



7* -2 

27^22 



7* -2 

2^33 



7, 



(36) 
(37) 



All of these points except for the Kasner circle lie in 11(1). 
5.1 Summary of stability properties 

We summarise the stability of those equilibrium points for the case 7^ > 2 (ie 
the anisotropic fluid is ultra-stiff). This was the situation in which the isotropic 
Friedmann-Lemaitre model was shown to be the single attractor on approach 
to the initial singularity in past studies [3], [6]. The particular cases in which 
at least one of the inequalities in p4p or ([55)1 on a,r or become equalities 
will be considered later. Note that our time derivative, r, is defined so that 
the universe is getting larger towards the future. Therefore, in the context of 
contracting universe, we are interested in the past asymptotic behaviour with 
respect to r. In terms of stability, the unstable equilibrium points are important. 

5.1.1 a < 0, r > 1 

In this case both Ai and A2 are located outside physical domain. J- is unstable 
and JC has stable arc and an arc of saddle points. 
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5.1.2 a<0, r<l 

A2 is not physical. T is unstable and K, is stable. Ai is unstable in the invariant 
set f2 = but a saddle point in the interior of the entire space. 

5.1.3 0<a<l, a<r^<l 

All equilibrium points lie in the physical domain. T is unstable in the invariant 
set Bianchi I and a saddle point in general. Ai is another saddle point. /C is 
stable. A2 turns out to be unstable and will be identified as the past attractor 
of the entire system. 

5.1.4 < a < 1, < a 

A2 is outside physical domain. Instead Ai becomes unstable and replaces the 
role of A2 in the previous case. 

5.1.5 0<a<l, r>l 

In this case Ai is outside the domain while A2 is physical and unstable. The 
difference is a part of Kasner circle becomes saddle. 

5.1.6 a > 1, r < 1 

The isotropic fluid is no longer stiff and A2 becomes non-physical. is the 
global future attractor, fC is saddle and Ai is the past attractor. 

5.1.7 a > 1, r > 1 

The situation is the same as the previous one except for non-physical Ai. A 
part of Kasner circle becomes the past attractor. 

5.2 Monotone functions 

Finally, we list some monotone functions crucial to understand the asymptotic 
behaviour. 

(/) Monotone for r 9^ in B^I)/ (n(I) U B(I) U {S = 0}). 

Increasing for < (j) — < n and decreasing for — 7r<(/> — 0<O. 

fl Monotone decreasing for a < in Ba{1) /A{T). 

For a = 0, it is scmi-monotonc decreasing and 11' = iff S = 0. 

From the monotonicity of (j) we conclude that it is sufficient to look at 11(1) 
and B{T) in order to determine the asymptotic behaviour . In particular, for 
a < l. r < 1, the past attractor (if it exists) must lie in 11(1) because the Kasner 
circle is stable. 
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5.3 Invariant set A(l) 

For completeness and to get a flavour of the dynamics involving anisotropic 
fluid, we consider the invariant set A(l). Setting = 0, the equations read 

S'± = {q-2)Y.±+3V±Z, (38) 

q = 2(S^ + E2_)+ 1(3^^_2)Z, (39) 

z = i-(i;^ + s2_). (40) 

We can easily eliminate q and Z to obtain 

'1 



(41) 



E± = -3(1 - - I 
To see the asymptotic behavior, we recast (|4T|) into 

We can see immediately that S± — is monotone increasing or decreasing 
according as 7^ > 2 or 7^ < 2. The equilibrium point Ai is a past attractor for 
7, > 2 and a future attractor for 7^ < 2. 

For 7^ = 2, we can find the exact solution. This is a critical case in which 
the equilibrium point disappears. From (|39p and (|40l) . we have q = 2. The 
dynamical equations then read 

E'± = 3V±Z, 
Z' = -6{r+T.+ +r^Y.-)Z, 

and give 

-^(P+I]+ + V-^-f = -{Vl + Vl). 
From this equation, we derive 

iV+T.+ +V-^-f = -{'Pl+'P-)Z + A^, (42) 

where A^ is a positive constant. Using this first integral, we arrive at the exact 
solution 



V-i-A 

^± = Vllvl tanh3Ar±y^^i3 (43) 



42 1 

Z = 5 (44) 

+ cosh^ 3At ^ ' 

where B is another integration constant. There is only one independent param- 
eter because equation (|40p serves as a constraint 

{Vl+V^)B^ = 1. (45) 



vi + vi 
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The general solution for 7 = 2 corresponds to 

a cx Ji+^(^33-7'22)]r ^^^^^.^;4V^3\ (46) 

b cx e[i+*(^"-^^^)]" cosh^?i+^l+^3\ 3At (47) 

c cx e [^+* (^^^ ^] " cosh +^1+^3^3 3At. (48) 

This solution is past and future asymptotic to Kasner vacuum (connecting two 
different Kasner points). Kasner exponents are given by 

iP33-V,,)± . (49) 

' 11 + ' 22 "T ' 33 





1 


B 


V\ = 


3^ 


3^3 




1 


B 


P2 = 


3^ 


3\/3 




1 


B 


P3 = 


3^ 


3^3 



SAV: 



{v.2~v,,)± (51) 

' 11 + '22 + '33 

where plus sign is for future and minus for past. 
5.4 Invariant set 11(1) 

On this subset, the governing equations are reduced to 

(52) 

(53) 
(54) 

(55) 

(56) 

They exhibit significant similarity to the set of equations for Bianchi II. In fact, 
the mathematical structure of the equations is the same for both systems and 
we can employ the techniques developed for Bianchi II by Collins [IS] to analyze 
our system. 

First, by using the constraints we eliminate q and Z: 

n' = 3(7, - 2) [a(l - f7) - E] f7 = X, (57) 
E' = ^(7,-2)[a(l-17)-E2 + l-a]E-^(7,-2)r(l-17-E2X58) 

= y, (59) 

For a > 0, let us define the function / , which was introduced by Collins HB], 
by 

/ = n-i (i-n-T,^) \ (60) 





= (q- 


2)S-^(7*-2)rZ, 


n' 


= (2g 


~37 + 2)f^, 


z' 


= [2?- 


-37, + 2 + 3(7, -2)rE] 


q 


= 2E2 


+ i(37-2)f^+i(37,-2)Z, 


1 


= fl + 


Z + E^. 
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Then, we have 



^ifX) + ^ifY) ^ -a)>0 (61) 

in the interior of 11(1). Hence / is a Dulac function for the system and we con- 
clude that there is no periodic or recurrent orbit in 11(1) by Dulac's theorem. 
If a < 0, we can just take — / as Dulac's function and are lead to the same 
conclusion. 

Combined with our knowledge of the eigenvalues for the equilibrium points, 
we can now draw conclusions about the past asymptotic behaviour of the sys- 
tem. However, special attention is necessary for the marginal cases with zero 
eigenvalues (see ref [TT| for detailed discussion of this situation in the absence 
of ultra-critical fluids) and they will be considered in the next section. 

5.5 Some special cases 

5.5.1 r = l: 

The equilibrium point Ai is identified with the Kasner equilibrium point <j) = 6. 
To determine the stability of this point we need to carry out a second-order 
analysis. This does not affect the conclusion that H(I) is the past- asymptotic 
set and Dulac's theorem holds for it. Therefore, the past attractor is either A2 
or Ai for < a < 1, and it is for a < 0. 

5.5.2 a = 0: 

Again, H(I) is past asymptotic set. On H(I) we can see from equation l|57p 
that f2 is monotone increasing unless S = 0. Therefore F is the past attractor 
although linearisation around it produces a zero eigenvalue. It is also clear that 
Ai is a saddle point and A2 is identical to J- . 

5.5.3 a = : 

This specialisation ensures Ai — A2 and the situation is the same as in the 
general case. 

5.5.4 a = l: 

This is the case where 7 = 2. The entire Bil) is a set of equilibrium points with 
two zero-eigenvalues. There are two sub-cases. If r < 1, Ai is the past attractor 
while A2 is unphysical. For r > 1, Ai is unphysical and A2 is identical to one 
of the equilibrium points in B{1)\ some of them have positive eigenvalues and $7 
is monotone increasing. Thus all orbits connect two of these equilibrium points. 
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5.6 Summary 

For the past asymptotic behaviour of we have found the following: 

• For a < 0, is the past attractor. 

• For < a < min{r^, 1}, A2 is the past attractor. 

• For a > r^,r < 1, Ai is the past attractor. 

• For a > l,r > 1, the Kasner circle is the past asymptotic set. 

• For a = 1, r > 1, a part of Kasner disc is the past asymptotic set. 

We give a plot in the parameter space of a by (Figure [1]) indicating the 
past asymptotic behaviour and showing phase portraits of the solutions in the 
invariant set 11(1) (Figure [H for some interesting cases. 



Figure 1: The different attractors that are approached near the singularity for 
different choices of the parameters a and r when 7^ > 2. The whole parameter 
space is divided into four regions. The past attractor is indicated for each of 
them by the symbol defined in the text. The isotropic Friedmann attractor is 
T; the other attractors are anisotropic. 



« A 




1 




A 



I 
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Figure 2: Phase portraits for different values of r and a. The Friedniann solution 
F is at (0,0) and the Kasner equilibrium points are located on (±1,0). The 
critical points Ai and A2 are also indicated in the figures. 
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6 Discussion 



We have studied the behaviour of an anisotropic universe containing ultra-stiff 
fluids with isotropic and anisotropic pressures. We find that the addition of 
anisotropic ultra-stiff pressures with principal pressures that can exceed p com- 
pletely changes the results obtained when only an ultra-stiff perfect ifuid with 
isotropic pressiire (p > p) is present. Most notably, the isotropic Fricidmann 
universe is no longer the stable early-time attractor solution as the initial sin- 
gularity is approached and the effects of the anisotropic pressures lead to an 
anisotropic Kasncr-like expansion near the singiilarity. We woiild expect that 
more general Bianchi type universes with ultra-stiff anisotropic pressures would 
lead to further types of anisotropic attractor but the Friedmann singularity 
would remain unstable. Any attempt to model the evolution of physical quan- 
tities through a singular or non-singular bounce in a cyclic cosmology will need 
to be re-evaluated in the light of these results. 

Acknowledement: K. Yamamoto acknowledges support from the Cam- 
bridge International Scholarship Scheme. 
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